## [1] "Census Tract" "Total Population"
## [3] "California County" "ZIP"
## [5] "Approximate Location" "Longitude"
## [7] "Latitude" "CES 4.0 Score"
## [9] "CES 4.0 Percentile" "CES 4.0 Percentile Range"
## [11] "Ozone" "Ozone Pctl"
## [13] "PM2.5" "PM2.5 Pctl"
## [15] "Diesel PM" "Diesel PM Pctl"
## [17] "Drinking Water" "Drinking Water Pctl"
## [19] "Lead" "Lead Pctl"
## [21] "Pesticides" "Pesticides Pctl"
## [23] "Tox. Release" "Tox. Release Pctl"
## [25] "Traffic" "Traffic Pctl"
## [27] "Cleanup Sites" "Cleanup Sites Pctl"
## [29] "Groundwater Threats" "Groundwater Threats Pctl"
## [31] "Haz. Waste" "Haz. Waste Pctl"
## [33] "Imp. Water Bodies" "Imp. Water Bodies Pctl"
## [35] "Solid Waste" "Solid Waste Pctl"
## [37] "Pollution Burden" "Pollution Burden Score"
## [39] "Pollution Burden Pctl" "Asthma"
## [41] "Asthma Pctl" "Low Birth Weight"
## [43] "Low Birth Weight Pctl" "Cardiovascular Disease"
## [45] "Cardiovascular Disease Pctl" "Education"
## [47] "Education Pctl" "Linguistic Isolation"
## [49] "Linguistic Isolation Pctl" "Poverty"
## [51] "Poverty Pctl" "Unemployment"
## [53] "Unemployment Pctl" "Housing Burden"
## [55] "Housing Burden Pctl" "Pop. Char."
## [57] "Pop. Char. Score" "Pop. Char. Pctl"
## [1] 0
## [1] 0
##
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|======= | 9%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|========= | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 15%
|
|=========== | 16%
|
|============ | 16%
|
|============ | 17%
|
|============ | 18%
|
|============= | 18%
|
|============= | 19%
|
|============== | 19%
|
|============== | 20%
|
|============== | 21%
|
|=============== | 21%
|
|=============== | 22%
|
|================ | 22%
|
|================ | 23%
|
|================ | 24%
|
|================= | 24%
|
|================= | 25%
|
|================== | 25%
|
|================== | 26%
|
|=================== | 26%
|
|=================== | 27%
|
|=================== | 28%
|
|==================== | 28%
|
|==================== | 29%
|
|===================== | 29%
|
|===================== | 30%
|
|===================== | 31%
|
|====================== | 31%
|
|====================== | 32%
|
|======================= | 32%
|
|======================= | 33%
|
|======================= | 34%
|
|======================== | 34%
|
|======================== | 35%
|
|========================= | 35%
|
|========================= | 36%
|
|========================== | 37%
|
|========================== | 38%
|
|=========================== | 38%
|
|=========================== | 39%
|
|============================ | 39%
|
|============================ | 40%
|
|============================ | 41%
|
|============================= | 41%
|
|============================= | 42%
|
|============================== | 42%
|
|============================== | 43%
|
|============================== | 44%
|
|=============================== | 44%
|
|================================ | 45%
|
|================================ | 46%
|
|================================= | 46%
|
|================================= | 47%
|
|================================= | 48%
|
|================================== | 48%
|
|================================== | 49%
|
|=================================== | 49%
|
|=================================== | 50%
|
|=================================== | 51%
|
|==================================== | 51%
|
|==================================== | 52%
|
|===================================== | 52%
|
|===================================== | 53%
|
|====================================== | 54%
|
|====================================== | 55%
|
|======================================= | 55%
|
|======================================= | 56%
|
|======================================== | 57%
|
|======================================== | 58%
|
|========================================= | 58%
|
|========================================= | 59%
|
|========================================== | 59%
|
|========================================== | 60%
|
|========================================== | 61%
|
|=========================================== | 61%
|
|=========================================== | 62%
|
|============================================ | 62%
|
|============================================ | 63%
|
|============================================= | 64%
|
|============================================= | 65%
|
|============================================== | 65%
|
|============================================== | 66%
|
|=============================================== | 66%
|
|=============================================== | 67%
|
|=============================================== | 68%
|
|================================================ | 68%
|
|================================================ | 69%
|
|================================================= | 69%
|
|================================================= | 70%
|
|================================================= | 71%
|
|================================================== | 71%
|
|================================================== | 72%
|
|=================================================== | 72%
|
|=================================================== | 73%
|
|=================================================== | 74%
|
|==================================================== | 74%
|
|==================================================== | 75%
|
|===================================================== | 75%
|
|===================================================== | 76%
|
|====================================================== | 76%
|
|====================================================== | 77%
|
|====================================================== | 78%
|
|======================================================= | 78%
|
|======================================================= | 79%
|
|======================================================== | 79%
|
|======================================================== | 80%
|
|======================================================== | 81%
|
|========================================================= | 81%
|
|========================================================= | 82%
|
|========================================================== | 83%
|
|=========================================================== | 84%
|
|=========================================================== | 85%
|
|============================================================ | 85%
|
|============================================================ | 86%
|
|============================================================= | 86%
|
|============================================================= | 87%
|
|============================================================= | 88%
|
|============================================================== | 88%
|
|============================================================== | 89%
|
|=============================================================== | 89%
|
|=============================================================== | 90%
|
|=============================================================== | 91%
|
|================================================================ | 91%
|
|================================================================= | 92%
|
|================================================================= | 93%
|
|================================================================= | 94%
|
|================================================================== | 94%
|
|================================================================== | 95%
|
|=================================================================== | 95%
|
|=================================================================== | 96%
|
|==================================================================== | 96%
|
|==================================================================== | 97%
|
|==================================================================== | 98%
|
|===================================================================== | 98%
|
|===================================================================== | 99%
|
|======================================================================| 99%
|
|======================================================================| 100%
Data is reduced to contain only Santa Clara data to avoid knitting problem. Map of PM2.5 shows the annual mean PM2.5 concentrations based on 2019 ACS population estimates in census tracts
Map of Asthma shows the age-adjusted rate of emergency department visits for asthma
The fitted line seems to fit most of the points, but on the upper half of the plot, there are points showing high asthma-caused emergency department visits that are not as related to PM2.5 concentration.
##
## Call:
## lm(formula = Asthma ~ PM2.5, data = ces4_map)
##
## Residuals:
## Min 1Q Median 3Q Max
## -27.346 -13.335 -5.208 9.570 156.830
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -51.086 15.353 -3.327 0.000964 ***
## PM2.5 10.218 1.863 5.484 7.71e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 19.62 on 370 degrees of freedom
## Multiple R-squared: 0.07517, Adjusted R-squared: 0.07267
## F-statistic: 30.08 on 1 and 370 DF, p-value: 7.714e-08
An increase of 1 unit in PM2.5 concentration is associated with an increase of 10.218 asthma-caused emergency room visit; 7.5% of the variation in Asthma is explained by the variation in PM2.5
The residual plot is significantly skewed to the left, maybe a curve would fit better than a straight line.
Try the log transformation to solve the skewness problem. the fitted line looks better than the previous one.
##
## Call:
## lm(formula = log(Asthma) ~ PM2.5, data = ces4_map)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.70323 -0.35769 0.01364 0.41749 1.74734
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.9074 0.4631 1.959 0.0508 .
## PM2.5 0.2931 0.0562 5.215 3.07e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5919 on 370 degrees of freedom
## Multiple R-squared: 0.06847, Adjusted R-squared: 0.06595
## F-statistic: 27.2 on 1 and 370 DF, p-value: 3.065e-07
Now we can interpret that with with 1% increase in PM2.5 concentration, the age-adjusted asthma-caused emergence room visit will increase by about 29%. 6.8% of the percent variation in Asthma is explained by the variation in PM2.5 Now the skewness problem is solved, it looks more like a normal distribution around 0.
Here the blue are the positive residuals and red are negative residuals. The most negative residuals is around Stanford University with a -1.7 residual. Negative residuals means that there is an over-estimation in Asthma estimation. My guess is that there are mainly students living in the area. Students in general are younger and healthier and potentially has no money to go to emergency room if not desperately needed to.